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Abstract. We apply optimal control theory to determine the shortest time in which 
an energy eigenstate of a weakly anharmonic oscillator can be created under the 
practical constraint of linear driving. We show that the optimal pulses are beatings 
of mostly the transition frequencies for the transitions up to the desired state and 
the next leakage level. The time of a shortest possible pulse for a given nonlinearity 
scale with the nonlinearity parameter (S as a power law 6~°' with a c:^ 0.73 ± 0.029. 
This is a qualitative improvement relative to the value a = 1 suggested by a simple 
Landau-Zener argument. 
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1. Introduction 



Physics is fundamentally quantum. On the other hand, in practice, many objects 
and processes appear to follow the laws of classical physics. This includes objects 
coupled to large baths with manydegrees of freedom, hence suffering from decoherence. 
Furthermore, in objects of large mass, quantum phenomena occur at an unobservably 
short de Broglie wavelength [1] . Harmonic oscillators, on the other hand, even if they are 
only weakly damped and have low mass and thus must be in the quantum regime, may 
appear close to classical: These linear systems are semiclassical, meaning that in their 
ground state and thermal equilibrium state, as well as in any states created from these 
by an arbitrary linear drive, correspond to the classical limit small quantum fluctuations, 
essentially following Ehrenfest's theorem. The smallness of these fluctuations and their 
competition with thermal fluctuations makes it hard to verify the quantum nature of 
these systems. A way out is to add an extra nonlinear ingredient to create a nonclassical 
state, such as a Fock state (an energy eigenstate of the oscillator) or a squeezed state 
(a state with quantum fluctuations in one variable below the standard Heisenberg 
limit). Examples include coupling the oscillator to a nonlinear object such as a qubit 
[21 El m El El [7j. In this work, we focus on a different approach, which is exploiting 
a small nonlinearity of the oscillator potential. These nonlinearities occur, e.g., in 
nanomechanical systems [HI E] , Josephson circuits [IDl IH [121 [El [E] , optics [151 [Ml [Ej , 
ion traps [TH [19], and molecules [20]. These systems that are rather large in size on 
an atomic scale are in particular important in this context, because this would enable 
to show the quantum character of center-of-mass variables of mesos- and macroscopic 
objects. Also, Fock states are useful for the demonstration of short-range quantum 
communication on a chip. 

The plan of the paper is as follows: we are first going to define our nonlinear 
oscillator Hamiltonian and outline the difficulty of creating Fock states in this system, 
as well as the GRAPE method used to find optimized ways to prepare these Fock states. 
We will then present optimized pulses, identifying a beating structure that can easily be 
connected to the energy spectrum. We are finally discussing the intrinsic speed limits 
for Fock state preparation. 



2. Setting and approach 

2.1. Hamiltonian 



H=^ + ^x^ + 5^x^ S>0 (1) 



We study a weakly nonlinear oscillator, the Duffing oscillator, with Hamiltonian 

2 2 2 2 

\ -x^ + 5 — 

2m 2 3^ 

where x and p are coordinate and momentum, respectively, where uoq is the resonance 
frequency in the limit 5 — > 0. Even though we focus on the hard case, 5 > 0, we 
expect similar conclusions to hold for the soft case. This Hamiltonian can be made 
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dimensionless with raising and lowering operators analogous to the harmonic oscillator 
[2T] . It then reads 

H = hiUo(^a^a + ^^ +^{a + a^)\ (2) 

By expanding the nonlinearity and ordering its terms, we can split the Hamiltonian as 
H = Ho{n) + V into a terms Hq that contains only the number operator h = a^a that 
will dominate perturbation theory for small 6/hu <C 1 and corresponds to the rotating 
wave approximation, and a term V that goes beyond that approximation. We find 

Hoin) = huJo(n + ^^ +^{2n^ + 2n + l). (3) 

V =^{a^^ + 6a^^ + 4a\a^ + a^^)a + 6a^ + a''). (4) 

Hq can be re-parameterized to a simpler form using u = llJq + 6, while dropping a 
constant energy shift 



Hn = n 



ujfi + -6n{n — 1] 



(5) 



where u is the frequency of the |0) ^ |1) transition. H can be diagonalized 
straightforwardly and its eigenstates \n), n = 0,1,2,... are close to the harmonic 
oscillator eigenstates, see fig. [H at small 6 and the eigenenergies are given by the 
polynomial deriving from eq. En — HqIu). We can note that for small 6, the energy 
splittings between adjacent levels are only slightly different from each other. Denoting 
the transition frequencies as cj-y = {Ei — Ej)/h, this means that uJn+i,n = uj + 6n only 
has a weak n dependence by 5 -C cu. 

In the most common physical situation, external driving is given by 

H^{t) = F{t)x = hf{t){a + a^) (6) 

corresponding to current drive in the Josephson laser field in the ion trap case 

and an external force in the nanomechanical case. Here, f{t) is proportional to the force 
F [21]. Given that the eigenstates are almost harmonic oscillator eigenstates, this drive 
predominantly, in the sense of an approximate selection rule, couples adjacent levels, n 
and n + 1. 

This setting outlines two difficulties in the task of preparing Fock states starting 
from the ground state using a straightforward application of He- Firstly, there is no 
direct transition from the ground state to any state higher than |1) unless 6 is large 
enough to make the eigenstate corrections due to V, eq. (jl]), significant. Secondly, even 
for preparing the first excited state, the approximate Liouvillian degeneracy (i.e., the 
weak 77,-dependence of ujn+i,n), prohibits short pulses: A short resonant vr-pulse of length 
Tr has, by the energy-time uncertainty principle of Fourier transform, a bandwidth of 
ub — Tc/Tr. If ljb becomes comparable to the small parameter \ujn+i,n — ^n,n-i\ = 
a pulse applied to u!n,n-i will have a significant Fourier component at Un+i^n- Hence, 
driving the transition from n—1 to n will inevitably also drive n to n+1 inducing leakage 
to the next higher state n + 1. This could be overcome by a very long vr-pulse that keeps 
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Figure 1. Ratio of the eigenergies of the RWA Hamiltonian, eq. [5] and exact 
eigenergies of the anharmonic oscillator for the first nine energy levels, assuming 
nonlinearities 6 = 0.01 and 6 — 0.1. 



00 B small, however, in a realistic setting, this will compete with energy relaxation back 
to the ground state at a rate Ti, effectively limiting ub ^ l/^i- 

2.2. GRAPE 

Given this conundrum, we resort to optimal control to find a way to prepare higher 
Fock states. We use the GRAPE algorithm [22]. This is an optimal control algorithm 
that uses open-loop control for optimizing control fields. The algorithm is: given an 
initial state a target state {ipf) and a time period for the state transfer T, such 
that \ip(to)) = and \ip(to + T)) = \ipf). The search starts from an initial guess for the 
controls, which for us is always taken to be a constant-envelope resonant vr-pulse with 
frequency u, i.e., f{t) = Qr cosut. Qr turns out to be the Rabi frequency of the 1 
transition. The time T is sliced into intervals of length St = T/N with between 
101 and 401 and pulse amphtudes are assumed to be constant across the intervals. As 
a performance index, we use 

^=\{H^f)\'. (7) 
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If fi designates the value of the control field / in eq. (jSD in time step i, it can be shown 
that 



Here, pj is the projector on the state in time step j and we see the structure of eq. ([6]). 
Here Xj = \ipj){ipj\ is the projector on the back-propagated target state 



Here, Ui is the propagator for the Schrodinger equation across time step j. The gradient 
eq. ([H]) can be very efficiently calculated this way and used for a gradient search for 
a locally optimal set of values fi describing the desired pulse shape. We temporarily 
increased the gradient by a factor of two for each time the pulse was found to be in an 
undesired local optimum. 

3. Numerical results 

We have analyzed this setting in a number of cases. In all cases, the initial state 
was chosen to be the ground state, \ipi) = |0) and the first, second, and third excited 
state where chosen as final states, \ipf) = \n = 1,2,3). We have truncated the energy 
spectrum at n = 10 levels and verified, that changing the number of levels did not lead 
to a discernible change in the pulse. In fig. [2] we show optimized pulses that display 
clear beating behavior. This beating becomes more complex as we go to higher states. 
More clearly, this structure can be analyzed by its Fourier transform, as shown in fig. [3l 
We observe a cluster of close discrete lines, confirming the observation that the pulses 
are beatings. The predominant frequencies in the |0) —>■ |2) transition are ooio = 1.12, 
the slightly higher uj2i = 1.24 and 00^2 = 1-36, with additional sidebands. Analyzing 
more transitions, we see the same picture, where a transition from |0) to \n) contains 
frequencies u!m+i,m for < m < n. 

This structure can be understood as follows: Due to the oscillator selection rule, 
the ladder of states needs to be climbed sequentially, it is not possible to go from \n) 
to \n + k) for any A; > 1, explaining the frequencies Uiq through uJn,n-i ^01 a transition 
^ n that correspond to driving vr-pulses on those respective transitions. Due to the 
bandwith issue discussed in section 12. H this also drives the transition to level n + 1, 
it is also necessary to drive the leakage transition leading out of |n), i.e., at ujn^n+i- A 
similar effect is seen in quantum gates for three-state systems |23l EH ESJ 

The state dynamics during the pulse can be seen in Fig. HI This analysis confirms 
the picture of sequential occupation of Fock states, climbing an energy ladder, with 
some intermediate leakage to higher state that corrects itself in the end. Next to these 
long-scale Rabi dynamics, there is a fast modulation on the scale of the average driving 
frequency typical for strong driving where the counter-rotating component of the drive 
that is not part of standard Rabi physics [21] becomes important. In the end of the 
pulse, populations of |1) and |3) are brought back to |2). 




(8) 



\^^) = U]^,---U, 



(9) 
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Figure 2. The amplitudes of three pulse sequences optimized by GRAPE for 
preparation of |1) (top), |2) (middle), and |3) (bottom) are shown. Each system had a 
Hilbert space consisting of 10 eigenstates, with S = 0.12. Each optimized pulse leads 
to a fidelity of 0.9999. 

4. Discussion 

As mentioned in tlie introduction, nonlinearity is a resource for the creation of Fock 
states. At very short pulse durations tg state transfer pulses must scale with a constant 
area, i.e., their amplitudes will scale as l/tg. Thus, if the pulse is too short the driving 
amplitude will be so high that the S-term in the Hamiltonian can be neglected and Fock 
state preparation will not be possible. We have numerically investigated the gate fidelity 
as a function of pulse duration and found the error to be roughly exponentially growing 
at short times, see fig. [5l We now investigate the minimal time to reach a fidelity of 
0.99999% for pulses |0) |1) and |1) |2), fig. [61 We see that this minimal time is 
a power law of the nonlinearity parameter 6; tmin oc 6°' with aoi = —0.73 ± 0.029 and 
= -0.90 ±0.031 

As a reference, we can construct the corresponding power law for non-optimized, 
single frequency Rabi pulses, starting from the |0) — >■ |1) transfer. Transforming to 
the frame rotating with the driving frequency, the effect of the pulse on the unwanted 
|1) —>■ |2) transition is to switch a transition matrix element on and off. This matrix 
element has to be compared to the detuning between the transition \lu2i — cuio I ~ (5. In 
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Figure 3. Fourier specrum for preparing |2}. To the left and riglit of the zero frequency 
is the imaginary and real coefficients, respectively, of the decomposition. The main 
peaks in each pulse spectrum is labeled. The transition frequencies between |0) and 
|1) {uj = 1.12) are apparent in each pulse. The transition frequency 0^21 — 1-24 and 
'-^32 = 1-36 is apparent in the pulse preparing the |2) and |3) states. 



order to create a Rabi pulse of fixed area vr, the maximal amplitude of the transition 
matrix elements scales as 1/tg, and its rate of change thus scales as In order to not 
occupy the |2)-state in the end of the pulse, the dynamics of the |1) ^ |2)-transition 
must be fully adiabatic [21], i.e., following a standard result from Landau-Zener theory 
[26l [27l EH] , we demand -C 5^ leading to a minimal gate time of cx 1/S. We would 
not expect this scaling law to change for the — > 2 transition, which is ultimately a 
sequence of two transitions of length l/tg, only a different prefactor. 

Thus, we can conclude that our optimal pulses qualitatively extend the limits of 
Fock state preparation by changing the minimal time to a softer power law, from 1/5 
to i/(50.73±o.o29 Qj^g j^g^j^ conclude that this is due to quantum interference in the higher 
levels. The difference is made possible by temporarily occupying higher states and 
unpopulating them in the course of the pulse. It needs to be remarked, that our pulses 
quahtatively differ from quantum gate pulses, i.e., rotations of the full basis instead of 
changing state [23l [23] , where a 1/5 scaling of the minimal time is found. This suggests 
that state preparation, which is less constrained, can use quantum interference more 
efficiently. A related paper has been posted recently [29] that looks at preparation of 
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Figure 4. Populations of the different states during the pulse taking |0) into |2}. It is 
recognized, that |0) is first excited into |I} and that this one then goes into |2} whereas 
higher states are somewhat occupied during the pulse and then go to zero. Additional 
fast modulation is due to counter- rotating terms. 

|1) exclusively and finds different envelopes, probably due to additional constraints. 
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